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Abstract 

Under the Basel II standards, the Operational Risk (OpRisk) advanced measurement approach is not 

prescriptive regarding the class of statistical model utilised to undertake capital estimation. It has 

however become well accepted to utlise a Loss Distributional Approach (LDA) paradigm to model the 

individual OpRisk loss process corresponding to the Basel II Business line/event type. In this paper we 

derive a novel class of doubly stochastic a-stable family LDA models. These models provide the ability 

to capture the heavy tailed loss process typical of OpRisk whilst also providing analytic expressions for 

the compound process annual loss density and distributions as well as the aggregated compound process 

annual loss models. In particular we develop models of the annual loss process in two scenarios. The first 

scenario considers the loss process with a stochastic intensity parameter, resulting in an inhomogeneous 

compound Poisson processes annually. The resulting arrival process of losses under such a model will 

have independent counts over increments within the year. The second scenario considers discretization 

of the annual loss process into monthly increments with dependent time increments as captured by a 

Binomial process with a stochastic probability of success changing annually. Each of these models will 

be coupled under an LDA framework with heavy-tailed severity models comprised of a-stable severities 

for the loss amounts per loss event. In this paper we will derive anal3ftic results for the annual loss 

distribution density and distribution under each of these models and study their properties. 

Key words: Operational Risk, Loss Distributional Approach, Doubly stochastic Poisson Process, 
a-Stable, Basel II, Solvency II. 
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1. Motivation 



The modelling of operational risk (OpRisk) has taken a prominent place in financial quantitative mea- 
surement, this has occurred as a result of Basel II / Basel III regulatory requirements. As a result 
Operational risk (OpRisk) has become increasingly important to the banking industry to address these 
regulatory standards in light of globalization, complex financial products and changes in information 
technology, combined with a growing number of high-profile operational loss events worldwide. 

There was no widely accepted definition of OpRisk when the Basel Committee on Banking Supervision 
(BCBS) began discussions on OpRisk management at the end of the 1990s; see BOBS (1998). Often, 
OpRisk was defined as any risk not categorised as market or credit risk. Some banks defined it as the risk 
of loss arising from various types of human or technical error. Some earlier definitions can be found in a 
1997 survey conducted by the British Bankers Association (BBA). In January 2001, the Basel Committee 
on Banking Supervision issued a proposal for a New Basel Capital Accord (referred to as Basel II) where 
OpRisk was formally defined as a new category of risk, in addition to market and credit risks, attracting 
a capital charge. In the working paper BCBS (September 2001) on the regulatory treatment of OpRisk 
and in the revised Basel II framework BCBS (2004), the following definition of OpRisk was adopted. 
'^Operational risk is defined as the risk of loss resulting from inadequate or failed internal processes, 
people and systems or from external events. This definition includes legal risk but excludes strategic and 
reputational risk. " This definition did not change in the latest version of Basel II framework, BCBS 
(2006, p. 144). The International Actuarial Association, lAA (2004), has adopted the same definition of 
operational risk in the capital requirements for insurance companies. 

So OpRisk is indeed a broad category. The BCBS gives a further classification into seven event types of 
OpRisk (BCBS, 2006, Annex 9): Internal Fraud; External Fraud; Employment Practices and Workplace 
Safety; Clients, Products and Business Practices; Damage to Physical Assets; Business Disruption and 
System Failure; Execution, Delivery and Process Management; which serves to further illustrate the 
disparate nature of events in this class. Reputational and strategic risk do not fall under the OpRisk 
umbrella, and market and credit risks are treated separately, but almost any other event that may result 
in a loss to a bank, including legal action, may be termed OpRisk. 

Basel II considers three pillars, which, by their very nature, emphasize the importance of assessing, mod- 
elling and understanding OpRisk profiles. These three pillars are minimum capital requirements (refining 
and enhancing risk modelling frameworks), supervisory review of an institution's capital adequacy and 
internal assessment processes and market discipline, which deals with disclosure of information. Since 
this time, the discipline of OpRisk and its quantification have grown in prominence in the financial sector. 

To illustrate just how significant OpRisk can be to a financial institution, one only needs to consider 
OpRisk related events such as the 1995 Barings Bank loss of around 1.3 billion GBP; the 2001 Enron 
loss of around 2.2 billion USD; the 2004 National Australia Bank loss of 360m AUD; and the more 
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recent Societe Generale loss of 4.9 billion Euro. Each of which demonstrates the severity of loss process 
that should be modelled by OpRisk statistical models, providing strong motivation for heavy-tailed loss 
process models such as those involving infinite mean and variace as captured by the family of a-stable 
models considered in this paper. 

The impact that such significant losses have had on the financial industry and its perceived stability 
combined with the Basel II regulatory requirements have significantly changed the view, that financial 
institutions have regarding OpRisk. Under the three pillars of the Basel II agreement, set out in the 
framework, internationally active banks are required to set aside capital reserves against risk, to imple- 
ment risk management frameworks and processes for their continual review, and to adhere to certain 
disclosure requirements. 

Whilst many OpRisk events occur frequently and with low impact (indeed, are 'expected losses'), others 
are rare, and their impact may be as extreme as the total collapse of the bank. The modelling and 
development of methodology to capture, classify and understand properties of operational losses is a new 
research area in the banking and finance sector. 

There are three broad approaches, that a bank may use to calculate its minimal capital reserve, as 
specified in the first pillar of the Basel II agreement. They are known as Basic Indicator Approach, 
Alternative Standardized Approach and Advanced Measurement Approach (AM A). In this paper the 
approach considered is the AMA. AMA is of interest since it is the most advanced framework with 
regards to statistical modelling. 

A bank adopting the AMA must develop a comprehensive internal risk quantification system. This 
approach is the most fiexible from a quantitative perspective, as banks may use a variety of methods and 
models, they believe are most suitable for their operating environment and culture, provided they can 
convince the local regulator, [BCBS, 2006, pl50-2]. The key quantitative criteria is that a bank's models 
must sufficiently account for potentially high- impact rare events. In this paper we consider the idea of 
the loss distribution approach (LDA) which involves modelling the severity and frequency distributions 
over a predetermined time horizon, typically annual as specified in the APS115 section on soundness 
standards. 

The fitting of frequency and severity distributions, as opposed to simply fitting a single parametric annual 
loss distribution, involves making the mathematical choice of working with compound distributions. 
This would seem to complicate the matter, since it is well known, that for most situations, analytical 
expressions for the distribution of a compound random variable are not attainable. However, the special 
classes of a-Stable models developed in this paper overcome this complication. 

Typically, the reason for modelling severity and frequency distributions separately then constructing a 
compound process because some factors affect the frequency and others may affect the severity, see|13|. 
Some of the key points relating to why this is important in most practical settings are that the expected 
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number of operational losses will change as the company grows. Typically growth needs to be accounted 
for in forecasting the number of OpRisk losses in future years, based on previous years. This can easily 
be understood, when modelling is performed for frequency and severity separately. Economic inflationary 
effects can be directly factored into size of losses through scaling of the severity distribution. Insurance 
and the impacts of altering policy limits and excesses are easily understood by directly altering severity 
distributions. Changing recording thresholds for loss events and the impact this will have on the number 
of losses required to be recorded is transparent. 

The most popular choices for frequency distributions are Poisson, binomial and negative binomial. The 
typical choices of severity distribution include exponential, WeibuU, lognormal, generalized Pareto, the 
g-a.nd-h family of distributions [B], [2] and recently the a-Stable family [20] . 

The most important processes to model accurately are those, which have relatively infrequent losses. 
However, when these losses do occur they are distributed as a very heavy-tailed severity distribution. 
In particular we focus our analysis on the scenarios involving heavy-tailed severity models in the rare- 
event extreme consequence context. Thereby providing analysis of the loss process most likely to have 
significant consequences on a financial institution, those which may lead to ruin. This involves introducing 
to OpRisk modeling an important family of severity models, utilized in insurance claims reserving in [I], 
given by the a-stable severity model. This family of severity model are flexible enough to incorporate 
light-tailed Gaussian loss models through to infinite mean, infinite variance severity loss models such as 
the Cauchy model. 

There are many approaches, which can be used to fit and incorporate expert opinion / scenario analysis 
for these parametric distributions and the approach adopted by a bank will depend on the data source 
being modelled and how much confidence one has in the data source. After the best-fitting models are 
selected, these are combined to produce a compound process for the annual loss distribution. From this 
compound process, VaR and capital estimates may be derived. 

Once compound processes have been fitted for each business unit and risk type, the next step in the 
process is to aggregate these annual loss random variables for each individual {business line-event type} 
combination, and thus to obtain the institution-wide annual loss distribution. This paper will not ad- 
dress the issues associated with correlation and dependence modelling. For more information on typical 
approaches to introducing correlation in an aggregation process, including copula methods, correlation 
of frequency, severity or annual losses, see [S] and |16j . In the next section we present the details of the 
LDA modelling framework adopted in this paper. 

2. Loss Distributional Approach Model Specifications 

OpRisk LDA models are discussed widely in the literature; see e.g. [5], [4] [9], [22]. Under the LDA 
Basel II requirements, banks should quantify distributions for frequency and severity of OpRisk for each 
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business line and event type over a one-year time horizon. In this section we begin by presenting a generic 
LDA framework before presenting specific models that we will develop under the compound process a- 
Stable family. The standard LDA Basel II structures, involving an annual loss in a risk cell (business 
line/event type) modeled as a compound random variable, 

(2.1) 

s=l 

Here i = l,2,...,r, T+lin our framework is discrete time (in annual units) with T + 1 corresponding 
to the next year. The upper script j is used to identify the risk cell. The annual number of events N^: is 
a random variable distributed according to a frequency counting distribution P^-' -'(•), typically Poisson. 
The severities in year t are represented by random variables xi'^\t), s > 1, distributed according to a 
severity distribution P^-'^(-), typically lognormal. Severities represent actual loss amounts per event. 

To reflect both the nature of OpRisk data with extreme but rare events, the severity models selected 
for this analysis are chosen to exhibit extreme heavy-tails, with particular interest in distributions with 
infinite mean/ variance. To accomplish this we consider the family of a-Stable severity distributions. 
These particular models have been proposed as suitable models for insurance claims modeling and finance 
in previous papers, such as [7], [H], [S], [H] and [T7]- 

The study of the distribution of the annual loss process is one of the classical problems in risk theory. 
Closed-form solutions are typically not available for the distributions used in OpRisk. However with 
modern computer processing power, these distributions can be calculated virtually exactly using numerical 
algorithms, see approaches in |TS]. The easiest to implement is the Monte Carlo method. However, 
because it is typically slow, Panjer recursion and Fourier inversion techniques are widely used. Both 
have a long history, but their applications to computing very high quantiles of the compound distribution 
functions with high frequencies and heavy tails are only recent developments and various pitfalls exist. In 
this paper the class of OpRisk models developed will be proven to admit closed form analytic solutions 
to the density and distribution functions of the annual loss and of the total bank's loss, for independent 
risks, in year t which is calculated as 

Z, = J2zi'\ (2.2) 

where formally for OpRisk under the Basel II requirements J = 56 (seven event types times eight business 
lines). However, this may differ depending on the financial institution and type of problem. In this paper 
we will drop the index j for the risk cell unless required. 

The approach we consider in this paper utilizes the representation of the annual loss process for a given 
business line and event type combination given by a convolution. In particular, it is well known that 
the density and distribution of the sum of two independent continous random variables Yi ^ ^i(') and 
I2 ^ -^2(0 with densities /i(-) and /2(-) respectively, can be calculated via convolution as 

/yi+i'2(2/) = (A */2) (y) = / f2iy - yi)fi{yi)dyi 
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and 

Py.+yM = iFi*F2) {y) = J F2{y ~ yi)h{yi)dy, 

respectively. Here the notation /i * /2 denotes the convohition of /i and /2 functions as above and the 
notation Y ^ F{y) means a random variable has a distribution function F{y). Thus the distribution of 



the aggregated loss in Equation 2.1 can be calculated via convolutions as 

oo 

F{z) ^ Pr[Z <z\^^Pr[Z < z\N = n] Pr [N = n] 

(2.3) 

Here, = Pr[Xi + X2 + ■ ■ ■ + X „ < z] is the n-th convolution of J^(') calculated recursively as 

F^'^>{z)= ( i^("-i)*(z)(z-a;)/(x)da; 
Jo 

with F(")*(z) = 1 if z > and zero otherwise. 

Next, we present a special family of statistical models for the severity distribution which are well known 
to be closed under convolution, meaning the analytic solution to the recursive convolution integrals above 
are known in closed analytic parameteric form for both the distribution and densities. 



2.1. a-Stable Severity Distribution Models 

Considered as generalizations of the Gaussian distribution, a-Stable models are defined as the class 
of location-scale distributions which are closed under convolutions. As in [20] we restrict to the class of 
truncated a-stable models to ensure we only work with non-negative loss processes. In an OpRisk context, 
a-stable distributions possess several useful properties, including infinite mean and infinite variance, 
skewness and heavy tails and pT| . 

We assume the i-th loss of the j-th risk process in year t m a risk cell is a random variable with a-stable 
distribution, denoted by xl-'\t) ^ Sa (x; /?, 7, (S, 0). Where, Sa (x; /3, 7, (5, 0) denotes the univariate four 
parameter stable distribution family under parameterization S'(O), see algorithm in the Appendix and 
details contained in [TT] . 

The univariate a-stable distribution we consider is specified by four parameters: a e (0, 2] determining the 
rate of tail decay; /3 E [^li 1] determining the degree and sign of asymmetry (skewness); 7 > the scale 
(under some parameterizations); and (5 G M the location. The parameter a is termed the characteristic 
exponent, with small and large a implying heavy and light tails respectively. Gaussian (a — 2, j3 — 0), 
Cauchy (a = = 0) and Levy (a = 0.5, /3 — 1) distributions provide the only analytically tractable 
sub-members of this family. Except these special cases, in general a-stable models admit no closed- 
form expression for the density which can be evaluated point-wise, inference typically proceeds via the 
characteristic function, see discussions in ^17^. However, intractable to evaluate point-wise, importantly 
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for OpRisk applications, simulation of random variates is very efficient, see [3] and the algorithm provided 
in Appendix |8] 

From 12], a random variable X is said to have a stable distribution, 5(a, /?, 7, 5; 0), if its characteristic 
function has the following form: 

exp{-7"|6l|"(l + i/3(sign(6i)) tan(^)(|76l|i-" - 1)) + iSO} ii a ^ 1 
exp{-7|0|(l + z/?(f )(sign(0))ln(7|0|)) + tSO} if a = 1. 



E[exp{i9X)] = 



In the following Lemmas we present some fundamental basic facts about a-Stable random variables that 
will be required to establish the novel analytic closed form expressions for the annual loss LDA models 
we develop in this paper. These will be used to construct analytic exact binomial, negative binomial, 
Poisson and doubly stochastic Poisson mixture representations of the annual loss process for a bank under 
a-stable severity models with the required positive support. This will be achieved by considering a special 
sub-family of a-stable models. 

In addition we will provide an analytic expression for the tail distribution of these models and the 
properties of the ES in special cases of this LDA model. Finally we will utilise known asymptotic results 
to comment on the institution wide capital aggregation results as a result of models containing dominating 
risk processes such as those contained in the a-Stable family. These results will extend those developed 
for OpRisk insurance models recently in |20j . 

Lemma 1 If Y ^ S{a, /3,^,S;0), then for any a 7^ 0,5 G 3?, the transformation Z = aY + b is a 
scaled version of the a-stable distribution. That is Z ^ S{a, {sign{a)l3, \a\"f,ad + b;0). In addition, the 
characteristic functions, densities and distribution functions are jointly continuous in all four parameters 
(a, /3, 7, S) and in x. These results follow from \21}j and Proposition 1.16. 

Lemma 2 If for all i G {1, . . . , iV} one has random variables Xi ^ S{a, Pi, 7^, Si, 0) then the distribution 
of the linear combination, given N, is 

N 

Z = ^X, ~^(a,^,7,?;0) 

p ^ T^^.Paf J ^ f S.=i 5. + tan ^ p7 " Eti ft'T.) " ^ 1 ^^'^^ 

Eti 7f \ Eti + l log 7 - Pjlj log 1) a = 1 

This result follows from ff^l}/ Section 1.2, Property 1.2.1) and (JT^, Proposition 1.17). 

The property in Lemma 2, of closure under convolution for random variables with identical a parameter 
is lost under truncation to positive support. The implications of this are that in OpRisk the convolution 
property will in general be lost for general members of the truncated severity models {X > 0). This is 
with one notable exception given by the sub-family of Levy distributions. Examples of loss distributions 
for this sub-family are illustrated in Figure [T] as a function of the scale parameter 7. 

Lemma 3 If X ^ 5(0.5, 1, 7, (5; 0) this model specifies the sub-family of a-stable models with positive 



1=1 

N 
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real support x € [5,oo). The density and distribution functions are analytic and given respectively, for 
5 < X < oo, by 



fx{x) 



ix-6) 



3/2 



exp 



7 



2{x~S) 



, Fx{x) = erfc 



7 



2 {x - 6) 



The median is given by ji — ^ {erfc ^ (0.5)) and the mode M = where erfc{x) = 1 — erf{x) 
1 /o^ ^ * This result follows from JT^ (Chapter 1. p.5). 




X(t) - Loss 



Figure 1: Example of the Levy severity model as a function of scale parameter, with location 5 = 0. 



2.1.1. Estimation of a- stable severity models in OpRisk 

From the perspective of application of such a-stable models in OpRisk LDA settings, there are efficient 
parameter estimation routines available given OpRisk data; see a discussion on Bayesian approaches and 
review in |17] . There are also available numerical maximum likelihood based approaches, see |21| and 
the software package implementing such numerical procedures as discussed in 12 . In this section we 
discuss the practical perspective of utilising a Levy stable model for the severity distribution in an OpRisk 
context and how this would be justified in practical data analysis. Firstly, in an OpRisk LDA setting the 
severity distribution will only take positive support, hence the restriction of /3 = 1 will provide this result 
by ensuring a positive skew that results in support only on the positive real line. We now discuss the 
appropriateness in practice of considering the Levy class of a-stable models, which satisfy the setting in 
which 13 — 1 and importantly admit closed form expressions for the density and distribution functions to 
be which we demonstrate in this paper results in closed form expressions for the annual loss distribution. 

We illustrate this reasoning on the following example, where we consider N = 100 losses Xi.^ obtained 
from a a-stable model Sa=o.7 (1, 7 ~ 0.01, S = 0). We fit general a-stable models for an increasing number 
of annual losses n G {5, 10, 15, 100}, the estimated MLE parameters and associated 95% C.I. are provided 
in Table |2.1.1[ These results were obtained utilising the software avialable on John Nolans a-stable 
website 

\protect\vrule widthOpt\protect\href {http : / / acadeinic2 . ainerican . edu/\string~ jpnolan/ stable/ stable . html}{l 
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MLE parameter 
estimates (95%C./[-]) 


N = 5 


N = 10 


N = 15 


N = 100 


OL 


1.01 [0,2] 


0.41 [0,2] 


0.45 [0,2] 


0.60 [0.46,0.74] 




0.04 [-1,1] 


0.99 [-1,1] 


1.00 [-1,1] 


1.00 [1.0,1.0] 


7 


0.51 [-0.14,1.16] 


0.14 [-0.23, 1.02] 


0.12 [-0.18, 0.84] 


0.01 [0.005,0.015] 


5 


-0.34 [0.31,0.99] 


0.12 [-0.26, 0.83] 


0.13 [-0.12, 0.56] 


0.02 [0.016, 0.024] 



Table 1: MLE parameter estimates of Sa (/3,7,5) and 95% C.I. 



The largest number of observations considered in this example with N = 100 is already practically a 
very large number of observations in the context of OpRisk in which heavy-tailed severity models are 
suitable and we see in this example that even with 100 observed losses, fitting the heavy tailed a-stable 
model estimates accurately the parameters /3, 7 and d. However, estimation of the tail index, in this 
class of model represented by a, with precision is a well known challenging problem with many classes 
of heavy tailed severity model. We demonstrate this here, where the tail index 95% C.I remains wide 
even with A'' = 100 observations. Therefore, when working with heavy-taile models in the a-stable family 
for an OpRisk setting it is sensible to consider selecting the sub-class of a-stable models which admit a 
parameteric form, this occurs for the Levy sub-family of a-stable models as presented in Lemma 3. We 
will now consider the LDA framework under this sub-class of severity model and demonstrate in such 
classes we can obtain closed form expressions for the annual loss distribution in several settings. 

Next wc develop novel LDA models utilising the sub-family of a-Stable Levy family severity models. 
Wc begin by considering discrete time loss processes in which the total number of lossess is assessed in 
fixed intervals of time, such as qiiarterly or monthly time increments. In such cases wc may wish to 
consider non-independent counts per month and so can consider the annual loss distribution modelled by 
a frcnquency distribution with binomial distribution or with negative binomial distribution, depending 
on whether wc believe the variance in the number of losses exceeds the mean number of losses or is less 
than the mean number of losses. In the continous arrival process setting, we consider annual losses to 
arrive accroding to independent increments of time according to an exponential arrival process, resulting 
in a Poisson frequency distribution. Clearly under a Poisson based model we arc implicitly assuming the 
mean number of losses per year is equal to the variance in the number of annual losses, when the rate 
parameter represents the average number of annual losses. We will then extend this Poisson family of 
models to the doubly stochastic family of frequency distributions and derive closed form expressions for 
the annual loss distribution in such cases. 
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2.2. binomial and negative binomial a-Stable LDA Analytic Models 



In this section we consider the binomial and negative binomial models where we assume that the random 
variable for the total number of events per year are first given by standard binomial and negative binomial 
Processes, before extending these results to doubly stochastic processes. For the binomial proces we obtain 
the following result in Theorem 1. 

Theorem 1 The distribution of the annual loss process Z represented by a compound process model with 
LDA structure in which the frequency is N{t) ~ Bi{M,p) and the severity model 

Xi{t) ~ S'(0.5, 1,7,(5;0), then the exact density of the annual loss process can be expressed analytically as 
a mixture density comprised of a-stable components with binomial mixing weights for Nt > 0, 



M 



fz{z) = ^c:fipni-pr- 



\ 3/2 



exp 



In 



X I 



5n < Z < OO 



(2.5) 



with 

n 



0.5 



= 1 5n = ^5t + tan - - ^ 1, = + tan - (n^|7| - 717) 



and Fz (0) = Pr{Nt = 0) = {1 — p)'^ for N = 0. The exact form of the annual loss cumulative distribution 
function is also expressible in closed-form, 

/ 



M 



Pr {Z<z)= Fz{z) = ^ {pT (1 - P)"""" erfc 



n=l 



7n 



2 (2 - 5n) 



X I 



5n < Z < QO 



(2.6) 



This result follows directly from application of Lemma 1, Lemma 2 and Lemma 3. 

For the doubly stochastic binomial-Beta proces we can state the following result in Theorem 2. 

Theorem 2 The distribution of the annual loss process Z represented by a compound process model 
with LDA structure in which the frequency is N{t) Bi{M,p), p Be {a, (3) and the severity model 
Xi{t) 5(0.5, 1,7, 5; 0), then the exact density of the annual loss process can be expressed analytically as 
a mixture density comprised of a-stable components with Beta mixing weights for Nt > 0, 

M 1 

fz{z) = V / {pf (1 -p)^-" ^7^,P"-\^-pr-'dp 



M „l 



[z - 6r^ 



3/2 



exp 



B{a,py 

In 



[z - (5„) 



X I 



5n < Z < CO 



M 



= E 

n=l 



M\B{a + n,l3 + M ~ n) 
{M -n)\n\B{a,fi) 



z — 5, 



3/2 



exp 



X I 



5n < Z < (X> 



(2.7) 
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with 



'0.5 



^-n|7|°■^ /3„ = 1 (5„ = ^5i+tan^ 7„-^7j I =n(5 + tan^ (n2|7| -717) 

1 i=i V j=i 



and Fz{0) = Pr{Nt = 0) = ^B(a^^"* f^f N = 0. Note we also denote the standard Beta function by 
B{x,y) = jQt'^~^{l —tY^^dt. The exact form of the annual loss cumulative distribution function is also 
expressible in closed-form, 

I 



M 



try . i:, / X ^ M\B{a + n,l3 + M -n) , 
Pr {Z <z) = Fz{z) = 2^ — — erfc 



n=l 



{M - n)\n\B{a, p) 



In 



2 - (5„) 



X I 



6n < Z < 00 



(2.8) 



+ 



B{a,p + M) 



B{a,p) 

Proof The result for the convolution between the stable random variables, conditional on n, follows 
from application of Lemma 1, Lemma 2 and Lemma 3. The mixing weight of the compound process for 
the binomial-Beta distribution is derived by considering Bayes Theorem and the conjugacy property of 
the binomial-Beta model from which we know p{p\n) = Be{p; a + n, ^ + M — n), hence 



p{n\M,p)pip) Bi{n;M,p)Be{p;a,P) (M-n)!n! P"(l " P)'^""^ B(a j) ""' 



Pr(n) 



p{p\n) 



Be{p; a + n,P + M -n) 



p(Q + „_l)(j^_p-)0 + M-n-l) 

B(a+n,fi+M-n) 



(2.9) 



M\B{a + n,P + M -n) 
{M -n)\n\B{a,p) ' 



□ 



For the negative binomial proces we obtain the following result in Theorem 3. 

Theorem 3 The distribution of the annual loss process Z represented by a compound process model with 
LDA structure in which the frequency is N{t) ~ NB{r,p) and the severity model 

Xi{t) ^ 5(0.5, 1,7, 5; 0), then the exact density of the annual loss process can be expressed analytically as 
a mixture density comprised of a-stable components with negative binomial mixing weights for Nt > 0, 



/z(^) = ^Cr'^-i(l-p)(n)(p)^ 



with 



In 



'0.5 



n=l 



1 



{z - 5„) 



3/2 



exp 



(yZ - Sn^ 



X I 



6n < Z < 00 



(2.10) 



/3„ = 1 (5„ = +tanM 7„ - ^7j =nd + t&nj (n^|7| - 717) 



i=l 



and Fz{0) = Pr{Nt = 0) = (1 — p)" for N = 0. The exact form of the annual loss cumulative distribution 
function is also expressible in closed-form, 



Pr {Z<z) = Fz{z) = J2 C:^'^-' (1 - pY {pT erfc 

n=l 

+(i-pr- 



In 



^ 2 - 5n) 



5n < Z < CO 



(2.11) 
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For the doubly stochastic negative binomial-Beta proces we can state the following result in Theorem 4. 

Theorem 4 The distribution of the annual loss process Z represented by a compound process model 
with LDA structure in which the frequency is N{t) ^ NB{r,p), p Be{a,/3) and the severity model 
Xi(t) ^ 5(0.5, 1,7,(5;0), then the exact density of the annual loss process can be expressed analytically as 
a mixture density comprised of a-stable components with Beta mixing weights for Nt > 0, 

X I < z < ooj 



exp 



7n 



= E 



(n + r- + rn, /3 + n) 



(r - l)!n!p('»+'-"-i)B(a,/3) 



3/2 



exp 



X I < 2 < CXjj 



(2.12) 



with 



'0.5 



TT 

tan - I 7„ 



^ 7j ^nS + tan - (71^17! - nj) 



and Fz{0) — Pr{Nt = 0) = ^ ^{a-i) for N = 0. The exact form of the annual loss cumulative distribution 
function is also expressible in closed-form, 



Pr(Z < 2) = Fz(2) = }^ ^ m ^er/c 



n = l 



(r- l)!n!p(°=+™-i)B(a,/3) 



2 z - <5„ 



X I 5„ < z < 00 



(2.13) 



+ 



r,{"-l) 



Proof The result for the convolution between the stable random variables, conditional on n, follows 
from application of Lemma 1, Lemma 2 and Lemma 3. The mixing weight of the compound process 
for the negative binomial-Beta distribution is derived by considering Bayes Theorem and the conjugacy 
property of the binomial-Beta model from which we know p{p\n) — Be{p; a + rn, /3 + n), hence 



Pr(n) 



p{n\r,p)p{p) ^ NB{n-r,p)Be{p-a,P) ^ G^'' ' (1 " vT {pY " ° 5(1 j) 
p{p\n) Beip; a + rn, l3 + n) 

{n + r- l)!p("+'-i)B(a + rn,l3 + n) 



p(a + r„-l)(l_p)(g + , 

i5(a-t-rn,/3-t-n) 



(r- l)!n!p("+™-i)B(a,^) 



□ 



(2.14) 



2.3. Poisson and Doubly Stochastic Poisson-Gamma - a-Stable LDA Analytic Models 

In this section we begin by reviewing some results recently developed in [20] before extending these 
results to the class of doubly stochastic Poisson-Ganima-a-Stable LDA models. Hence, we begin with 
the standard Poisson process LDA model result given in Theorem 5. 
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Theorem 5 The distribution of the annual loss process Z represented by a compound process model with 
LDA structure in which the frequency is N{t) -Po(A) and the severity model 

Xi{t) ^ S{0.5, 1,7, 5; 0), then the exact density of the annual loss process can be expressed analytically as 
a mixture density comprised of a-stable components with Poisson mixing weights for > 0, 



fz{z) = J2^M-^) 



n=l 



A" 



n! 



X I 



6n < Z < OO 



(2.15) 



with 



In 



'0.5 



i=l 



i=i 



and Fz{Q) ~ Pr{Nt = 0) = cxp(— A) for N = 0. The exact form of the annual loss cumulative distribution 
function is also expressible in closed-form, 

/ 



Pr {Z<z)= Fz{z) = ^ exp(-A) — er/, 



n=l 



In 



2{z-5r 



X I 



5n < Z < OO 



+ exp(-A). (2.16) 



This result follows directly from application of Lemma 1, Lemma 2 and Lemma 3. 

We can now present results for the doubly stochastic Poisson-Gamma process LDA model, given in 
Theorem 6. 

Theorem 6 The distribution of the annual loss process Z represented by a doubly stochastic compound 
process model with LDA structure in which the frequency is N{t) ~ -Po(A) and A F (a, /3) and the 
severity model 

Xi{t) ~ S{0.5, 1,7, 5; 0), then the exact density of the annual loss process can be expressed analytically as 
a mixture density comprised of a-stable components with Poisson mixing weights for Nt > 0, 



exp(-A)^ ^ exp (-/?A) A^^dA 
n! 1 (a) 



fz{z)=f:\[ 

n=l 

(a-l)!n! \1 + 13 J \l + P 



2- {z-S.f' 



exp 



X I < « < ooj 



X I |^5„ < « < ooj 



(2.17) 



with 



-0.5 



1 Sn = ^Si+ta.nj I -'in ^Ij] =nS + tan j (n^|7| - n-y) 



and Fz{0) = Pr{Nt = 0) 



()3+l)° 
r(a) 



for N = 0. The exact form of the annual loss cumulative distribution 
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function is also expressible in closed-form, 



PTiZ<z) = Fziz) = J2 



{a + n- 1)! 




Sn < Z < CO 



(2.18) 

Next we present the proof of this result. 

Proof The resuh for the convolution between the stable random variables, conditional on n, follows from 
application of Lemma 1, Lemma 2 and Lemma 3. The mixing weight of the compound process for the 
Poisson-Gamma distribution is derived by considering Bayes Theorem and the conjugacy property of the 
Poisson-Gamma model from which we know p{X\n) — T{X; a + n, /3 + 1), hence 



Pr(n) 



p{n\X)p{X) Po{n;X)r{X;a,(3) 



p{X\n) 
r(a + n) 



r(A;a + n,/3 + l) x a+n-i^-fa+iu 

r(a+n) 



(2.19) 



/3" 



(a + n- 1)! / P 



T{a)n\ (1 + /3)" 



(a-l)!n! \l + P J \l + P 



Note, this mixing weight is then a negative binomial probability with r = a and p — as presented in 



p.53 of Ho] 



□ 



Next we will also study the truncation error in these settings for single loss processes. This will be 
followed by a consideration of some simulation studies comparing the properties of the capital estimation 
under the binomial, negative binomial and Poisson models and the doubly stochastic models developed. 



3. Approximation of Infinite Sum Representations by Truncation 

In this section wc consider analytic expressions for the approximation error of the infinite sum LDA model 
representations in Theorems 3,4,5,6. To achieve this we will build on the results developed in [20], which 
utilsed the properties in Lemma 4. 

Lemma 4 Given a random variable for the severity of a loss X ^ S{a, /?, 7, 6; 0) then as x ^ 00 one can 
write the limiting tail distribution 

P(X > x) - 7"c„(l + /3)2:"", as x^ 00 

, (3-1) 

fx{x\a,fi,^,5]Q) a^°'ca{l + I3)x ("+^', as x^ 00. 

where = sin (^) ^^^^ ■ This result follows from fflS^. Theorem 1.12) 

We may now use Lemma 4 to help determine a truncation index for the infinite sums in Theorems 3,4,5 
and 6 as show in Theorem 7. To achieve this we will first make an approximation in which we consider 
the index of the summation to be a continous variable, as demonstrated in for example [Section 3.1] of 
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[T5] . In addition, in this paper we will work with the analytic compound process expression for the tail 
asymptotic, and ensure that after truncation, the error in the finite expansion of the tail asympotic is 
bounded to a desired precision. More precisely, the asymptotic tail probability of the compound process 
will be established to have a desired approximation error under the truncation derived. This is not 
unique, and equally we could have chosen other distributional properties such as the median that was 
used in |2D]. 

Theorem 7 Given the LDA structure in which the frequency is distributed generically according to 
N{t) ^ /(•) and the severity model Xi{t) ^ S(0.5, 1, 7, S; 0), the exact expressions for the tail asymptotic 
of the resulting infinite mixture is given by, 

00 00 
Pr{Z{t) > zg) ^ 2a:-"-5co.5 ^ Pr{N{t) - n)j"„-^ = Ct„; ^ W„ 

n—l n—1 

where Zq is an upper tail quantile of the annual loss distribution and it will be assumed for simplicity that 
In^'^ = X]r=i ^t'^ ~ n|7|'^'^. These results allows us to determine a unique maximum of ^ \og(Wn) — 
corresponding to the term n with the maximum contribution to the LDA compound processes tail probabil- 
ity for each model. Note this can be determined trivially for each model by searching for a sign change in 
the function for each integer value of n and so is very efficient to solve since only requires integer search. 
The specified equations for each model are given as follows: 
Poisson-a-Stable. 

Under this model we have that Wn = exp(~A) i^,^__iy which results in the function ^log(M^„) = 

producing log(A) — log(ri) ^ 5^ = 0. 

Doubly Stochastic Poisson-Gamma-a-Stable. 

Under this model we have that Wn ~ ItT''^ ^'o^")!??!' ( T+/3 ) ( ) ^hich results in the function log( W„) 
producing log(n + a - I) + ^^^^^ + ^^^^ - log(n) - 1 - 5^ - log(/3 + 1) = 0. 
negative binomial-a-Stable. 

Under this model we have that Wn = {n+r-i-ny.{ny. ~ P)'' (p)" results in the function 

£ log(VK„) = producing log(n + r - 1) + ^^^^ + - 1 - log(n) log(p) = 0. 

Doubly Stochastic negative binomial- Beta-a- Stable. 



Under this model we have that Wn = ^"^(r-V)rri!p(°+^"^i"B(™'/f)"'"""' ""^^^"^^ results in the function 

^ log(W^„) = producing, 

, , , {n + r + 0.5) , , , r(a + rn- 0.5) , ,^ , (;9 + n - 0.5) 

log 71 + r) + ^- + r log a + rn) + r — ^ + log /? + n) + ^ 

1M . o N (r + l)(a + rn + /3 + n-0.5) , , , (?i + 1 + 0.5) , , , „ 

- r + 1 log a + rn + /? + n - ^ - log n + 1 - , - r logp = 0. 

(a + rn + p + n) (n+1) 

(3.2) 

We can now find Wl < Wq < Wjj such that the truncated sum approximation is sufficiently accurate for 
use in the expression of the annual loss distribution of the compound process. Hence, we can bound the 
number of terms in the mixture representations to not change by a precision amount e. This is expressed 
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analytically as 

^ ,3/2 exp I ^ ~ \ I (3-3) 

with Nl and Njj are determined from max(VFL < e^"^^Wo, 1) and Wjj < e^'^^Wg. The derivation of 
these results are found in Appendix [9j We will demonstrate simulated examples for these results on the 
truncation in the following results and discussion section. 



fz{z) = ^ Pr{Nt=n) 



4. Simulation Studies 

In this section we undertake several simulation studies, for all studies we will utilise the parameter settings 
in Table|4] In addition, the calculated values for N]^ and Njj are provided for each of these models, utlising 
the results derived in Theorem 7 in Section [3j For each of the frequency distributions, we considered two 
settings, low frequency and high frequency. In addition, all simulations were completed in Matlab with 
vectorised code on an SGI cluster in which the compute node utilised was an Altix XE320 with two Intel 
Xeon X5472 (quad core S.OGHz) CPU's, the source code available at: 

\protect\vrule widthOpt\protect\href {http : //web .maths .unsw. edu. au/\string~peterga/ index_f iles/QRSLab/Tei 



LDA Model 


Frequency Parameters 




Severity Parameters 


Truncation 




A 


n 


P 


r 


7 


6 


Nl 


Nu 


Low Frequency Example 


binomial-Levy 




12 


0.1 




0.01 








M 


binomial-Beta-Levy 




12 


Be(l,5) 




0.01 








M 


negative binomial-Levy 






0.1 


2 


0.01 








19 


negative binomial-Beta-Levy 






Be(l,5) 


2 


0.01 








20 


Poisson-Levy 


0.1 








0.01 








11 


Poisson-Gamma-Levy 


Ga(l,0.1) 








0.01 








171 


High Frequency Example 


binomial-Levy 




12 


0.6 




0.01 








M 


binomial-Beta-Levy 




12 


Bc(5,l) 




0.01 








M 


negative binomial-Lcvy 






0.6 


10 


0.01 








120 


negative binomial-Beta-Levy 






Be(5,l) 


10 


0.01 








169 


Poisson-Levy 


10 








0.01 








49 


Poisson-Gamma-Levy 


Ga(l,10) 








0.01 








171 



Table 2: Parameter settings for the models utilised in the simulations. 
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4-.1. Analysis of Model Estimation under Exact Truncated Solutions versus Monte Carlo 

In this section we first illustrate that the derived annual loss compound process distributions in Theorems 
1 to Theorem 6 are correct. This is achieved by undertaking a simulation study in which we simulate from 
each of the models developed first via a standard Monte Carlo procedure, utilising T = 1, 000, 000 annual 
simulated years for each LDA model. We record the simulation time as well as the estimated Monte Carlo 
error obtained by blocking the simulated years into sets of 50,000 samples, used to calculated the standard 
error in estimates obtained from the simulated annual losses. We then estimate the empirical distribution 
function and plot this along with the estimatd standard error from the Monte Carlo simulation. Next we 
evaluate each of the models using truncations of the exact expressions derived in Theorem 1 through to 
Theorem 6, where we utilise = 1 and Njj — 1, 000. The evaluations are done at 200 points from Z — 1 
to Z — 200 which for the parameter settings in Table |4] covered sufficient range for the evaluation to 
include quantiles representing regulatory capital VaR's at 99.5%. We then plot the evaluated distribution 
functions versus the Monte Carlo estimated empirical distributions functions for each of the compound 
process models. 

In Figure [2] we present the results of this simtulation for the Annual loss compound process with Levy 
severity distribution and the binomial and the doubly stochastic binomial-Beta models derived in Theorem 
1 and Theorem 2. The distribution functions are presented for both the low and high frequency examples. 

In Figure [3] we present the results of this simtulation for the Annual loss compound process with Levy 
severity distribution and the negative binomial and the doubly stochastic negative binomial-Beta models 
derived in Theorem 5 and Theorem 6. The distribution funcstions are presented for both the low and 
high frequency examples. 

In Figure [4] we present the results of this simtulation for the Annual loss compound process with Levy 
severity distribution and the Poisson and the doubly stochastic Poisson-Gamma models derived in The- 
orem 5 and Theorem 6. The distribution funcstions are presented for both the low and high frequency 
examples. 

5. Analysis of Truncation Error and Computational Saving versus Monte Carlo 

Monte Carlo techniques are widely used in the generation of compound processes and especially in the 
estimation of the annual loss distribution and the associated risk measures of Value at Risk (VaR) and 
Expected Shortfall (ES). In this section we fist study the truncation error when applied to the models 
derived in Theorem 1 through to Theorem 6 and then we demonstrate the accuracy of the derived models 
under the estimated truncations Nl and Nu . Finally, we compared the simulation times for each method 
when we very conservatively set = 1 and Nu — 1,000, well in excess of the required estimates 
truncations for these examples as shown below. 
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-Monte Carlo Empirical CDF 
Monte Cario +- 1 std. error 
-Ecact Solution (N^= 1, = 1000) 



— Monte Carlo Empirical CDF 
Monte Carlo +- 1 std. error 
... Exact Solution (N =1,N =1000) 



-Monte Carlo Empirical CDF 
■ Monte Carlo +- 1 std. error 
-Exact Solution (N =1,N = 1000) 



— Monte Carlo Empirical CDF 
Monte Carlo +- 1 std. error 
- . - Exact Solution (N^^ = 1 , = 1 000) 



Figure 2: Comparison of cpirical CDF from T = 1,000,000 simulated annual years (dashed lines), versus Truncated 
Nu = 1000 Exact Expressions (solid line). Top Left - Low frequency Poisson-Levy model Top Right - High frequency 
Poisson-Levy model Bottom Left - Low frequency Poisson-Gamma-Levy model Bottom Right - High frequency Poisson- 
Gamma-Levy model 

We demonstrate the estimated CDF and MSE for different levels of truncation Nu e {1, 2, ... , 1,000} 
in Figm'e[5]for the binomial and doubly stochastic binomial-Beta examples. Figure |6] shows the negative 
binomial and doubly stochastic negative binomial-Beta examples and Figure [7] shows the Poisson and 
doubly stochastis Poisson-Gamma examples. 

These results demonstrate that in all cases, the mean squared error falls rapidly for small numbers of 
upper truncation index Njj. This confirms the utility of utilising these truncations, which are extremely 
efficient simulation wise, where in Table [3] we see that even when using an excessive number of terms 
Njj — 1,000 we make very significant computational savings and accuracy improvements over the Monte 
Carlo equivalent simulations. 

Next we present a general result which provides a closed form distributional result for the combined loss 
process for multiple risk processes, each under one of the LDA model structures presented. 

6. Multiple Risk Process Aggregation 

In this section we demonstrate how to obtain the distribution of the combined loss processes under each of 
the models developed by presenting a general closed form result for the aggregation of two risk processes 
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-Monte Carlo Empirioal CDF 
" Monte Carlo +- 1 std. error 
-Exaot Solution (N =1,N = 1000) 



-Monte Carlo Empirical CDF 
Monte Carlo +- 1 std. error 
. Exact Solution (N =1,N =1000) 



-Monte Carlo Empiricai CDF 
■ Monte Carlo +- 1 std. error 
-Exact Solution (N =1,N = 1000) 



-Monte Carlo Empirical CDF 
Monte Carlo +- 1 std. error 
-Exact Solution (N|_= 1, N^= 1000) 



Figure 3: Comparison of cpirical CDF from T = 1,000,000 simulated annual years (dashed lines), versus Truncated 
Nu = 1000 Exact Expressions (solid line). Top Left - Low frequency Poisson-Levy model Top Right - High frequency 
Poisson-Levy model Bottom Left - Low frequency Poisson-Gamma-Levy model Bottom Right - High frequency Poisson- 
Gamma-Levy model 

which is trivially extended to any number of risk processes. Then we illustrate this result with an example 
using two Poisson processes LDA models. 

Theorem 8 The distribution of the annual loss process represented by multiple risks, eg. j G {1,2}, in 
a LDA compound process structure in which the frequency is N^^'^(t) ^ F{9^-^^) and the severity model 
Xl'\t) - 5(0.5, 1,7(J'),(5(^);0), can be expressed analytically as a mixture distribution comprised of a- 
stable components with mixing weights specified by the chosen frequency distribution. Define the total 
annual loss random variable as, 

z,=z«+zp)= y: e ^fw- (6-1) 

Furthermore, if each of these loss processes has an insurance mitigation applied under an ILP policy 
with top cover limits, S^^"^ > TCL^^"^ and (5^^-* > TCL'^^\ then the closed-form analytic expression for the 
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-Monte Carlo Empirical CDF 
Monte Carlo +-1 std. error 
-Exact Solution (N|^ = 1,N^ = 1000) 



-Monte Carlo Empirical CDF 
■Monte Carlo +- 1 std. error 
-Exact Solution (N|^= 1, N^^ = 1000) 



-Monte Carlo Empirical CDF 
■Monte Carlo +-1 std. error 
-Exact Solution (N =1,N = 1000] 



-Monte Carlo Empirical CDF 
■Monte Carlo +-1 std. error 
- Exact Solution (N =1,N =1000) 



Figure 4: Comparison of cpirical CDF from T = 1,000,000 simulated annual years (dashed lines), versus Truncated 
Nij = 1000 Exact Expressions (solid line). Top Left - Low frequency Poisson-Levy model Top Right - High frequency 
Poisson-Levy model Bottom Left - Low frequency Poisson-Gamma-Levy model Bottom Right - High frequency Poisson- 
Gamma-Levy model 

annual loss density for N^^'^ -\- A^^ > is given by 



oo oo 



m—Q n—0 



2tt 



(z (^nm) 



3/2 



exp 



2 Z-Sr, 



X I 



^0.5 

/ nm 



IL 111 I IL lib 



(2) 

7fc 



i=l 



fc=l 



^Jd) + ^^(2) + tan ^ (^(7^17(1^-^ + m|7('T-^)' - ^7^'^ - ^7(2)^ , 



(6.2) 



and /zt(0) — Pr{N^^'' + N^"^^ — 0). The exact form of the insurance mitigated annual loss cumulative 
distribution function is also expressible in closed-form, 

C30 00 

Pr {Zt <z)= Fz^z) = E E ^'■(^t = m)Pr{N^^'^ = n)erfc 
+ Pr{N^^^ +N^^^ = 0). 



A 



I X I [z > 0] 

2{Z-Snm] I (6.3) 
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Annual Loss Distribution versus Truncatic 



Mean Square Error of Truncation versus "Truncation Index [SJ 




- "^u = 
rsij^ = 30 

- N, , = 1 OOO 




upper Truncation Index [s^^ 



Annual Loss Distribution versus Truncation is^ Mean Square Error of Truncation versus Truncation Index hjlj 




=15 

- = 

, = 30 

- N, , = 1 OOO 



1 



Upper Truncation Index Isj^ 



Figure 5: Comparison of high frequency Poisson-Levy LDA model epirical CDF versus Truncation Threshold Njj - Exact 
Expression {Nu = 1000) (solid line). Left - Estimated CDF versus Njj Right - Mean Square Error versus Njj 



6.1. Example 1: Poisson-a-Stable Bivariate Risk Aggregation 

In the case of two Poisson process LDA models one would obtain the result given by 
f g„p(-A...-A.^.)'^"')"'f)" 



m— n— 



/ nm 



^71- / ? V' 

n m 



exp 



2 z-d„ 



(2) 

,7fe 



\ 3 = 1 fe=l 



X I z > (5„ 



(6.4) 



22 





Figure 6: Comparison of high frequency Poisson-Levy LDA model epirical CDF versus Truncation Threshold Njj - Exact 
Expression {Nu = 1000) (solid line). Left - Estimated CDF versus Njj Right - Mean Square Error versus Njj 

and fzTi^) — Pr(7Vf'-^'' + N^'^'' = 0) = cxp(— Ai — A2). The exact form of the insurance mitigated annual 
loss cumulative distribution function is also expressible in closed-form, 

/ I ^ 

xl[z>0] 



Py{Zt <z)^ Fz^z) = y y exp(-A(i) - A^^)) ^ > } ^ > erfc 

m— n— 

-h exp(-Ai - A2). 



In 



2 z-S„ 



(6.5) 



7. Conclusions 



This paper provided novel closed form analytic results for the several interesting classes of OpRisk LDA 
models based on the family of doubly stochastic a-stable LDA models. Thereby providing models with 
the ability to capture the heavy tailed loss process typical of OpRisk whilst also providing analytic 
expressions for the compound process annual loss density and distributions as well as the aggregated 
compound process annual loss models. 

Models for the annual loss process were developed in two scenarios, the first based on a loss process 
treated as continous in time throughout the year with a stochastic intensity parameter, resulting in 
an inhomogeneous compound Poisson processes annually. The second scenario considered discretization 
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Figure 7: Comparison of high frequency Poisson-Levy LDA model epirical CDF versus Truncation Threshold Njj - Exact 
Expression {Nu = 1000) (solid line). Left - Estimated CDF versus Njj Right - Mean Square Error versus Njj 

of the annual loss process into monthly increments with dependent time increments as captured by a 
binomial process with a stochastic probability of success changing annually. 

The resulting LDA annual loss distribution exact representations derived were then analysed and com- 
pared to a standard Monte Carlo estimates obtained after simulating 1, 000, 000 annual years. The results 
demonstrated how accurate the exact expressions were in both low and high frequency scenarios even 
with a small number of terms in the expansions. In addition, we studied the mean square error of the 
estated distribution functions in each of the six models developed as a function of the truncation. 

Finally, we demonstrated a novel multiple risk result showing how these resuls could be generalized to 
the multiple risk process setting to allow one to obtain expressions for the annual loss process of an 
institutions based on individual loss processes from one of the six doubly stochastic families of models 
developed. 
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Model Exact 


Truncatx-d (Nt = l.Nt = 1.000) 


Alontc Carlo (T = 1.000.000) 




Low Frequency Examples 
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6.1 (sec) 


10 (min) 
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1.2 (sec) 


8 (min) 
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Poisson-Levy 


1.2 (sec) 


8 (min) 


Poisson-Gamma-Levy 
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7 (min) 
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1.2 (sec) 


5 (min) 


Poisson-Levy 


1.2 (sec) 
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1.2 (sec) 


6 (min) 



Table 3: Simulation times for exa^t solutions versus Monte Carlo estimates. 
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8. Appendix 1 



Simulation of a univariate random variable from a general a-stable model can be achieved based on the 
following steps. 

1. Sample W ~ Exp(l) 

2. Sample U - Uniform f] 

3. Aply transfromation to obtain sample y split into a composite: 



(cost/) a 

(b) y = I [(I + I3u) tanu - /31n 



cos(ii— a(n+Sa ./3 )) 



if a^^ 1 



if a = 1 



^ arctan (/3tan (^)). In this setting y will be 



where Sa^fs = (l + /J^tan^ (^)) ° and S„,/3 
distributed by Sa (/?, 1, 0) 
4. Apply the transformation to obtain sample y = jy + S with parameters Sa (/3, 7, 5) 



9. Appendix 2 

Proof of Theorem 7 Given the LDA structure in which the frequency is distributed generically ac- 
cording to N{t) /(•) and the severity model Xi{t) 5(0.5, 1, 7, d; 0), the exact expressions for the tail 
asymptotic of the resulting infinite mixture is given by, 

00 00 
Pr {Z{t) > Zg) - 2x-° 'co.5 ^ Pr{N{t) = n)^°-' = Cxaii ^ 

n=l n=l 

where Zg is an upper tail quantile of the annual loss distribution and it will be assumed for simplicity that 
—0.5 _ ^^^^ ^0.5 _ j.j|^|0.5 Tiiese results allows us to determine a unique maximum of ^ log(W„) = 
corresponding to the term n with the maximum contribution to the LDA compound processes tail 
probability for each model as folows: 
Poisson-g-Stable. 

Under this model we have N{t) -Po(A) which results in the term W„ = |7|°'^ exp(— A) ■ Now we 
approximate this function as being continous in n and we differentiate and equate to zero, ^ log(W„) = 0, 
to find the mode of the terms in the series being summed. We first apply Stirling's approximation to 
n\ = r(n + 1) w \/27rn"+°-5e-" to give 

logW^„ = 0.51og(|7|) - A + nlog(A) -log((n- 1)!) 
= 0.51og(|7|) - A + nlog(A) - log(r(n)) 

w 0.51og(|7|) - A + nlog(A) - (0.51og(27r) + (n + 0.5) log(n) - n) 
^ ^ [log Wn] = log(A) - log(n) - ^ 
^ n large [log W„] « log(A) - log(n). 



27 



Equating this to zero and solving for n, tells us that if A > 1, then the maximum occurs at approx n = X. 
If A < 1, then max occurs at approx n = 1. □ 

Doubly Stochastic Poisson-Gamma-a-Stable. 

Under this model we have N{t) -Po(A) and A F (a, P) which results in the term W„ = |7|°'^ ^'{,^-V)hi\ (t+/3 
Now we approximate this function as being continous in n and we differentiate and equate to zero, 
^ log(W„) = 0, to find the mode of the terms in the series being summed. Again we apply Stirling's 
approximation r(n + 1) « \/2TTn"~^^'^e~"' to give 

log Wn = 0.5 log(|7|) + log (r(a + n)) - log (r(a)) - log (r(n + 1)) + a log(/3) - (a + n) log(/3 + 1) 
w 0.5 log(|7|) + 0.5 log(27r) + {n + a- 0.5) log(n + a - 1) - a + 1 - log (r(a)) - 0.5 log(27r) 
- (n + 0.5) log(n) + a log(/3) - (a + n) log(/3 + 1) 

Equating this to zero and solving for n, ie. solving for the roots log(2) — log(n) — 5^ = can be done by 
an explicit search, which produces n = 2. □ 

negative binomial-a-Stable. 

Under this model we have N{t) ~ NB{r,p) which results in the term 

Wn = 171°'^ (^n+r-i-ny.{ny. (-*■ ~ pY (p)"- Now wc approximate this function as being continous in n and 
we differentiate and equate to zero, ^ log(W„) = 0, to find the mode of the terms in the series being 
summed. We first apply Stirling's approximation to n! = T{n + 1) « y/2Trn'^~^^'^e~" to give 

log W„ = 0.51og(|7|) + log((n + r- 1)!) - log((r - 1)!) - log((n)!) + rlog(l - p) + nlog(p) 
w 0.51og(|7|) + 0.51og(27r) + (n + r - 0.5) log(n + r - 1) - (n + r - 1) - log((r - 1)!) 
— 0.5 log(27r) — (n + 0.5) log(n) +n + r log(l — p) +n log(p) 

^ = + + n + r -0.5 + - ' - - ^ + 

Then for a given set of values r and p in the model we can equate to zero this expresion and solve for n. 

□ 

Doubly Stochastic negative binomial- Beta-a-Stable. 

Under this model we have N{t) NB{r,p) and p ~ Be{a,P) which results in the term Wn = 
|^|0.5 . We will now utilise the well known identity for the Beta function 

with respect to Gamma functions and a Stirling Approximation given by 

™a;-0.5yj/-0.5 



Now we approximate this function as being continous in n and we differentiate and equate to zero. 
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f- log(W^) = 0, to find the mode of the terms in the series being summed. 



dn 

\ogWn = 0.51og(|7|) +log(r(n + r)) + (a + r - l)\og{p) + \og{B{a + rn, (3 + n)) 

- log(r(r)) - log(r(n + 1)) - (a + rn - 1) log(39) - log(B(a, /?)) 

« 0.51og(|7|) + 0.51og(27r) + (n + r + 0.5) log(n + r) - n-r+{a + r- 1) log(p) 
+ 0.5 log(27r) + {a + rn- 0.5) log(a + rn) + (/3 + n - 0.5) log(/9 + n) 
-{a + rn + P + n- 0.5) log(a + rn + ^ + n) - log(r(r)) - 0.5 log(27r) 

- (n + 1 + 0.5) log(n + 1) + n + 1 - {a + rn - 1) log(p) - log{B{a, /3)) 

d r, 1 , / N (n + r + 0.5) , , , r(a + rn- 0.5) , s 

^ — log Wn] = log n + r) + ^- — -— ^ + r log a + rn + ^- — ^ ^ — + log^ + n 

an (n + r) (a + rn) 

(/3 + n-0.5) , , « N (r + l)(a + rn + /3 + n-0.5) 

iP + n) ^ ^ ^ ^ (a + rn + ^ + n) 

(n + 1 + 0.5) 

- log(n+ 1) (n + l) ''loSW 

Then for a given set of values r and p in the model we can equate to zero this expresion and solve for 
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